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Abstract. In Active Galactic Nuclei (AGNs) the presence of a star cluster around the central black hole can have 
several effects on the dynamics and the emission of the global system. In this paper we analyze the interaction of 
stellar atmospheres with a wind outflowing from the central region of the AGN nucleus. Even a small mass loss 
from stars, as well as possible star collisions, can give a non-negligible contribution in feeding matter into the AGN 
nuclear wind. Moreover, stellar mass loss can produce envelopes surrounding stars that turn out to be suitable for 
reproducing the observed emission from the Broad Line Region (BLR). In this framework, the envelope can be 
confined by the bow shock arising from the interaction between the expanding stellar atmosphere and the AGN 
nuclear wind. 
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1. Introduction 

It is generally accepted that the broad emission lines 
present in AGN spectra originate in numerous, small and 
relatively cold gas concentrations in the neighborhood of 
the central black hole, the so-called Broad Line Region 
(BLR). The observed line properties directly indicate the 
existence of such gas concentrations; however, the physical 
nature of these structures is still a matter of debate. Up 
to the present day, three possible classes of models have 
been envisaged: gas clouds, accretion disks, and stellar 
wind envelopes. Each scenario suffers from still unclari- 
f ied pr oblems and we refer to the recent review by Korista 
( |l999| ) for a presentation of the problem. 

In general, broad emission line regions in AGNs turn 
out to be characterized by a relatively narrow range in 
ionization parameter U, i.e., U = L/[47 rcr^n,, < hv > 
] ~ 0.0 1 - 1 (Alexander & Netzer |l994 Peterson |l997 , 
Netzer 1990 ). In this expression r is the distance from 
the central black hole, n* is a representative value of the 
number density of the line emitting gas, L is the ionizing 
radiation field luminosity from a central source, illuminat- 
ing the line emitting region, and < hv > is the average 
photon energy of the ionizing radiation field. A first esti- 
mate of the conditions of the line emitting gas, such as n* 
and U, can be obtained from an analysis of the emission 



spectrum; coupling this information with the evaluation of 
the luminosity and of the ionizing spectrum, the relation 
above, defining [/, can be inverted to give an estimate of 
the characteristic size of the BLR, or, better, the charac- 
teristic distance of the line emitting gas from the central 
continuum source, tblr (see Wandel et al. 1999). In a first, 
simplistic use of photoionization arguments, under the as- 
sumption that the ionizing continuum shape is quite sim- 
ilar and both and U are substantially the same for all 
the AGNs (see Kaspi et al. 2000), a relationship, between 
the characteristic BLR size and the AGN luminosity, of 
the type tblr oc L'^'^ would then be expected. Indeed, 
this would just be an order of magnitude evaluation, since 
there are reliable indications that the BLR material den- 
sity and the ionization parameter are not constant at all 
even across the BLR extension of a single AGN (see, e.g., 
Kaspi et al. 20001 ). Nonetheless, this was actually what re- 
verberation mapping studies of several AGNs seemed to 

1/2 

suggest until very r ecentl y, impl ying t blr — 0.1 L^g pc 
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(Netzer & Peterson |1997[ Kaspi |1997|) (where L^e is the 
AGN luminosity in units of 10'*^ erg/sec). In a reverbera- 
tion technique study of a sample of 17 quasars, combining 
their data with those available for Seyfert 1 galaxies, Kaspi 
et al. ( pOOOl ) propose a significantly different relationship 
between tblr and the luminosity of the AGN, namely 
''BLR L*^'^. Although the exponent does not differ much 
from the previous one, it is remarkable that the ensemble 
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of the data fit by the latter relation, are not fit by the one 
with an 0.5 exponent. 

In a previous paper (Pietrini & Torricelli 2000, here- 



after Paper I), motivated by recent observations support- 
ing radial outflows even in radio-quiet AGN's, such as 
Seyfert Is, (Crenshaw et al. 1999, Weymann et al. 1997), 



we have analysed the physical structure and characteris- 
tics of a global AGN outflow, presumably originating in 
the very central regions and expanding out to large dis- 
tances as a kind of background/connection for various ob- 
servational components of the AGN structure [ BLR, UV 
absorbers. X-ray "warm absorbers"]. A solution for the 
wind equations, accounting for the physical requirements 
typical of central regions of Seyfert-like AGN, turned out 
to be possible only under rather specific conditions. One of 
the resulting constraints is that the outflowing wind may 
exist only if its density is rather low. The only possibil- 
ity to increase the wind density is to feed matter into the 
wind itself in a certain range of distances from the nucleus 
(see Sect. 8.4 of Paper I for details). 

As far as the BLR cloud model is concerned, the pres- 
ence of a nuclear wind generally augments problems of 
cloud survival, rather than solving them, unless the clouds 
are somehow comoving with the outflow itself. In fact, one 
problem for the cloud existence is the drag force at work 
between the moving clouds and a non-comoving confin- 
ing medium. This drag force would rapidly disrupt the 
clouds. The observationally inferred information on the 
cloud kinematics indicates that this motion is not radially 
directed and it must be characterized by a more complex 
velocity field. Also, the possibility of a Keplerian pattern 
has been suggested by Wandel et al. (1999). Since, for 
the reasons explained above, our choice is to work within 
the scenario of a radially directed nuclear outflow, the 
cloud survival problem would inevitably be worsened by 
the presence of a radial nuclear wind. 

For the stellar wind envelopes, since they are continu- 
ously fed by the stellar wind, disruption is not a problem; 
another velocity component in their relative motion with 
respect to the interstellar medium does not constitute a 
problem either. On the contrary, we will show in this pa- 
per that this new component in the relative motion of the 
envelopes with respect to the surrounding medium can 
be a fundamental factor; in fact, it can contribute signif- 
icantly to confine the plasma envelopes via the formation 
of bow-shock fronts. 

The above considerations, together with the necessity 
of feeding matter in the AGN wind, induce us to invoke 
stellar envelopes as a possible interpretation of the nature 
of BLR emitting plasma. In addition, the virial assump- 
tion of Keplerian motions used to interpret the emission 



line width (Wandel et al. 1999) is in good agreement with 
the stellar origin of BLR emission. Indeed, this possibil- 
ity has been widely discussed in the literature. In par- 
ticular, after two pioneering works of Scoville & Norman 



(1988) and Kazanas (198£), the star model has been re- 
analyzed by Alexander & Netzer (1994, hereafter ANl), 
and Alexander & Netzer (1997, hereafter AN2), who in- 



troduce "bloated" stars as stars characterized by a very 
extended envelope, but, due to the specific AGN environ- 
ment, different from supergiant stars as known in the solar 
neighbourhood. We refer to ANl for a wider presentation 
of the scenario and its problems. 



In their papers, Alexander & Netzer (AN1,AN2) 
present different possibilities for confining "bloated" star 
envelopes and conclude that the more effective one is tidal 
disruption by the black hole. However, ANl and AN2 as- 
sume in their works that the stellar winds expand into 
vacuum, asserting that this does not constitute a limita- 
tion on the validity of the analysis. On the contrary, in 
this paper we analyze whether and how the presence of a 
nuclear AGN wind, such as the one described in Paper I, 
can have an influence on the BLR physical structure. In 
particular, by inserting mass losing stars in their specific 
environment and investigating the interaction of the stel- 
lar wind with the AGN wind, we introduce another con- 
fining mechanism for the star envelopes. We then compare 
it to both the tidal one and to the other possible confine- 
ment mechanisms investigated by ANl and AN2. Finally 
we analyze the consequences for BLR parameter values. 



In Sect. 2, we discuss confinement of the expanding 
stellar envelopes and introduce the bow-shock formation 
mechanism as a further means of defining an outer bound- 
ary for these envelopes. Section 3 is devoted to the con- 
nection and interaction between the AGN nuclear wind, 
as developed in Paper I, and the central star cluster, with 
colliding and mass losing stars, in terms of the resulting 
mass deposition into the AGN nuclear wind. In nSect. 4, 
we describe the physical parameters that turn out to be in- 
volved in the computation of the model, their significance 
in the outcome determination and the general properties 
of a stellar envelope necessary to produce the observed line 
emission. We then identify a sort of "observational test" 
for the results of a BLR model, through the comparison 
with typical values inferred from observations of global 
BLR parameters, such as the covering factor, the charac- 
teristic BLR radius, the ionization parameter, and the fact 
that no broad forbidden lines are observed. Section 5 dis- 
cusses how the identified general requirements put physi- 
cal limitations on our model parameters. The general re- 
sults of our model are presented in Sects. 5.1, 5.2, 5.3 for 
three distinct cases, that differ basically in the the choice 
of the mass loss rates for the various types of stars in the 
central cluster. Section 6 is devoted to general discussion 
of the main results, including a comparison with the works 
of other authors applying different models to infer the gen- 
eral properties of the BLR of a specific source, NGC 5548, 
which is well representative of the Seyfert 1 class of AGNs. 
Finally, in Sect. 7 we summarize our results and the gen- 
eral features that characterize our interpretation of the 
BLR, highlighting the qualities and the limitations of the 
model. 
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2. Confinement of the stellar envelopes: Bow 
shocks around moving stars versus other 
mechanisms 

When a mass-losing star moves rapidly through the in- 
terstellar medium, its wind interacts with the surround- 
ing plasma, giving rise to an elongated bow shock. The 
physics of these shock fronts has been widely studied in 
differ ent co ntexts (see e.g Perry & Dyson 1985| , Van Buren 
et al. 199C and references therein) for different purposes. 
In general, we expect that the expanding stellar wind can 
develop a shock front, in which the wind kinetic energy 
is partially converted and radiated away. Therefore, in a 
stationary configuration, the shock position is determined 
by the equilibrium of the external pressure and that of the 
stellar wind. The resulting shock geometry is a cometary- 
like tail elongated in the direction opposite to that of the 
star motion. 

In our case the situation is more complex since stars 
move around the AGN nucleus with a velocity V(r) and 
are embedded in a radial nuclear wind characterized by 
a number density n^{r) and expanding with velocity 
Vw(7') = v„ir)e,-. The AGN wind velocity Vv/(r), number 
density n^(r) and temperature T^{r) profiles are those 
computed through the resolution of the system of sta- 
tionary hydrodynamical equations, described in Paper I, 
including the appropriate mass deposition. The nuclear 
wind velocity undergoes the transition from sub-sonic 
to super-sonic at a distance rc from the AGN nucleus, 
where rc/r^ > l/4(c/cs)^, with Cg being the local sound 
speed, c the Hght speed, and = 2G'Mbh/c^ the black 
hole Schwarschild radius. 

For the stellar velocity we adopt a Keplerian motion 
(consistent with recent interpretations of emission line 
width (Wandel et al. |1999D): 



V{r) = 



GM{r) 



(1) 



where M{r) is the total mass (central black hole plus star 
cluster) within the radius r, assuming its distribution is 
spherically symmetric. The velocity of the stars is gener- 
ally supersonic inside the nuclear wind sonic radius, while, 
of course, is supersonic outside. Hence, in our case the 
bow shock will have a different geometry depending on 
the regime in which it is computed: comctary-like shocks 
with the tail in the direction opposite to the star motion 
(when the ram pressure contribution is induced by the 
star motion, mn'T-wV^^); or with the tail extending radi- 
ally outwards from the central black hole (when the ram 
pressure contribution is due to the AGN wind: mn^iw^J); 
or a combination of the two. It is beyond the scope of 
this paper to investigate the detailed physics and geom- 
etry of this configuration; instead, what is important in 
this context is the existence of a region around the star 
where the stellar wind can be confined by this mechanism. 
For the sake of simplicity we have chosen to represent this 
complex wind envelope region schematically as a quasi- 
spherical structure, so that its geometrical extension can 



be characterized simply by a radius defining the radial 
distance of the envelope boundary from the star. 

In order to derive this envelope radius, first we have 
to parameterize the physical properties of the expanding 
envelope, namely give a prescription for the density of the 
gas in the envelope itself and for its expansion velocity. 
Our simplifying choice is to adopt a radial power law be- 
haviour for the stellar atmosphere expansion velocity: 



R 

v^R) = f'o ( 



for R>R^ 



(2) 



where R, the distance from the star center, and i?*, the 
radius of the star, are both expressed in the same length 
units; this relation is taken as independent of the distance, 
r, of the star from the central black hole. Taking into 
account continuity equation for the envelope gas, we can 
derive its density as a function of R as well: 

■ , , -("+2) 

"-onv(P) 



R 



Once these quantities are defined, the order of mag- 
nitude of the distance from the star, i?bow, at which the 
bow shock can be formed, can be determined by using 
the standard methods, that is from the two equations de- 
scribing the stellar mass loss rate (i.e., continuity equation 
again) and the balance between the total pressure (ther- 
mal-f-ram pressure) of the stellar wind and that of the 
ambient medium, all evaluated at -Rbow^ 



A-h = 47r mn i?bow"-cnv(Pbow)w*(Pbow) 

P* + mHne„v(i?bow)wJ(Pbow) = -Pw + Pra 



(4) 
(5) 



Here we assume that P» = ncnvfcP*/^* and = 
n^kT^/fi^ are the thermal pressures of the star and of 
the AGN winds, respectively, where rii = Pi/mn (for 
i = env,w), with pi the mass density. We have taken 
P-* = Pw = 0.6, as expected for a fully ionized plasma 
with cosmic abundances (see Allen 1985| ); note that in the 
following we indicate both /i, and /^w as p,. Also, is the 
expansion velocity of the stellar atmosphere as given by 
Eq. (H) for R — i?bow The ram pressure of the external 
medium can be identified with the maximum of that due 
to the star motion and that due to the AGN wind: 

Pram = MoX [mH^w t/^ , TOnriw^w] = "iH?^wWMax• 

The two equations above (Eq. (^ and Eq. (||)) turn 
out to depend on the distance of the star (whose enve- 
lope we are describing) from the central black hole, r, 
since they contain physical quantities depending on r it- 
self. These are the ones defining the AGN nuclear wind 
{nyj{r),T^(r),v^{r)), and the Keplerian velocity of the 
star in the local gravitational field, as given by Eq. (|l|). 
Also, these quantities can be combined so as to define 
global envelope properties as functions of r as well. In 
fact, Eq. (0) can be rewritten as 



n*(r) 



/(Pb 



(6) 
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and, making use of Eq. (^, Eq. (||) becomes 



-Rbow('') 



(7) 



where indicates the value of the stellar atmosphere ve- 
locity at its boundary i?bow, as in Eq. (||). Note that in 
Eq. (|^) we have redefined ncnv(^bow), the density of the 
stellar envelope at its external boundary, as ex- 
plicitly showing its dependence on r, since its value does 
depend on the distance of the envelope from the black 
hole, as mentioned above. In the following, unless other- 
wise specified, we use n*(r) to indicate a reference number 
density value characterizing the stellar atmosphere. As for 
the stellar atmosphere gas temperature, here we assume 
for simplicity that it is constant, i.e. independent both of 
r and of R; this assumption and its validity within the 
present framework are briefly discussed in Sect. 5. The 
two relationships shown above define the density and the 
extension of the stellar wind envelope as functions both of 
the stellar wind velocity and temperature, and of the ex- 
ternal confining plasma parameters. Since the stellar wind 
velocity depends on i?bow , Eq- ^ is in general an implicit 
equation. 

As Alexander & Netzer (ANl, AN2) have shown, sev- 
eral other mechanisms in principle can be competing for 
the definition of the physical extension of the stellar wind 
envelope. In the following we briefly outline them, in or- 
der to compare the characteristic envelope radius values 
defined by each mechanism at any given radial distance r 
from the black hole with i?bow , *• e- the distance from the 
star at which the bow shock, due to the interaction of the 
stellar wind with the nuclear AGN outflow, would form. 
The various physical processes that are to be examined 
produce a stellar wind envelope extension that is charac- 
terized by a different dependence on the radial distance r 
from the central black hole. Therefore, for each value of r, 
we shall finally identify the effective extension of the enve- 
lope itself with the smallest radius obtained, thus selecting 
the most efficient mechanism for stellar wind confinement. 

A first effect we need to take into account is the pos- 
sibility of tidal disruption of the outer layers of the stellar 
envelope, due to the black hole gravitational pull. In this 
case the size of the stellar envelope is given by (see ANl) 



R 



tidal 



1/3 



(8) 



showing a simple linear dependence on the distance from 
the black hole. 

Another possible definition of an upper limit to the 
stellar wind extension can be determined by the reason- 
able assumption of the existence of a finite mass in the 
stellar wind. Again following ANl, we suppose that the 
total mass in the stellar wind must not exceed 0.2 M0; 
this leads to the condition 



ALdt = 



dR = 0.2Mr. 



(9) 



which, for a given stellar wind velocity field w*(i?), defines 
-Rmass as the maximum stellar wind size allowed for, due 
to the limit on the stellar wind mass content. Assuming 
that Eq. (|^) defines the velocity profile (i?) for the stellar 
wind, we can explicitly solve the integral to obtain 



(1 - a)0.2Moi;o 



i?* R 



1/(1-") 



(10) 



notice that in this specific case the stellar envelope exten- 
sion would not depend on the radial distance r from the 
central black hole. 

Comptonization effects, due to the central ionizing 
continuum heating the outer layers of the stellar envelope, 
should also be taken into account in the determination of 
the external boundary of the envelope itself. This process 
(see Kazanas 1989 , and ANl) induces complete ionization 
of the outer layers of the stellar wind, that become part 
of hot "Comptonized" phase. It is effective on the denser 
stellar wind down to layers in which the density of the 
envelope material reaches a critical value, ncomp, above 
which, that is deeper into the envelope, a "cool phase" 
is still allowed (Kazanas 1989| ; Krolik, McKee, & Tarter 
1981), and therefore the stellar wind survives as such. This 
density value is defined by a critical value, Sc, of the "ion- 
ization" parameter S = Lion/ (47rr^cnfcT) for the existence 
of a "cool" phase, for a given ionizing radiation field, and 
it is 



"-Comp(r) = — 



1 



1 



(11) 



This critical value of the density defines a boundary radius 
for the stellar wind, that we can call "Comptonization" ra- 
dius, -Rcomp, by setting ncnv{R) = 'T-Comp in the continuity 
equation (see Eq.(H)). We thus obtain 



^Co 



miiVoR:;"Lion 



(12) 



In equations (11) and (12) the quantities and S can- 
not be freely chosen, since they depend on each other, 
their specific relationship, r^(S), deriving from the anal- 
ysis of possible equilibrium configurations existing in a 
photoionized gas (Krolik, McKee, & Tarter 1981 ). Since 
T^CE.) depends on the assumed continuum radiation illu- 
minating the g well as on other heating mechanisms 
possibly present and on cooling processes, its resulting 
detailed shape turns out to be different in literature de- 
pending on the specific assumptions of the authors. Krolik 
( |1999D shows a T{E) curve, obtained with the composite 
continuum spectrum from Elvis et al. (1994), in which 
Sc — 10 and the corresponding cool phase limit tempera- 
ture is ~ 5 X 10^ K. In other studies (see Krofik 2002), the 



"cool" phase (T < lO^K) exists for ^ < 28 for t he contin- 
uum spectrum considered by Kaspi et al. (2001), whereas 
for Krolik & Kriss ( pOOl| ) it is present only for ^ < 12. 
The temperature corresponding to these upper limits for 
the ionization parameter S is around 5 x 10'*K for the two 
cases. It is somehow difficult to make a consistent choice 
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for the values of Sc and T* to be used in Eq. (12) without 
computing directly the specific equilibrium, therefore we 
choose to rely on the assumptions of one of the papers in 



characteristic number density n* can be derived from the 
continuity equation [Eq.(@) with i?ext the place of i?bow] 



literature, namely the one of Kazanas (1989). In this work, 
the author computes the critical density for the transition 
to the Comptonized phase (see Eq. (11)) taking the critical 
value of the ionization parameter as Sc = 10 and a tem- 
perature (T* in our notation) equal to 3 x 10^ K (Kazanas 
1989). Given the results of most recent studies and calcu- 



lations (some of which have been cited above), these values 
can be regarded as sort of "lower limits" for the possible 
range of critical ionization parameter and/or temperature. 
However, in the present context we maintain this choice 
as representative of the critical condition, since, owing to 
the functional dependence of the limiting radius, i?compj 
as expressed in Eq. (12), "lower limit" values for and 
the corresponding T^, guarantee a "lower limit" value for 
the Comptonization radius as well. In fact, we are evalu- 
ating this radius to compare the efficiency of this confine- 
ment mechanism with that of the other relevant ones. If 
indeed Comptonization turned out to be the effective con- 
finement mechanism for the stellar envelopes contributing 
to the line emission in the BLR, a more accurate anal- 
ysis of appropriate parameters to be used in the precise 
evaluation of -Rcomp would be required. Nevertheless, if, 
as it will be proved to be the case in all the conditions 
of interest for envelopes contributing to the construction 
of our BLR model, the physical extension of the stellar 
envelopes is not determined by -Rcomp, the present evalu- 
ation is fully representative. Also, it is important to note 
that, if this is the case, the Comptonization radius limita- 
tion itself is no longer physically relevant to the problem. 
In fact, outside the otherwise confined stellar envelope, 
the gas is mixed up with the surrounding nuclear wind 
and shares its physical properties. 

Finally, ANl mention another limit for the external 
radius of the stellar wind envelope, deriving from the re- 
quirement that the stellar wind velocity, when it is in- 
creasing outwards, does not reach values larger than 0.1c. 
This is a rather artificial constraint, defining a radius i?v, 
which, for our choice of the expansion velocity of the stel- 
lar atmosphere, is i?v — R* (O.Ic/wq)^^", and should apply 
only in extreme cases. 

In summary, we define the effective stellar envelope 
extension at each distance r from the central black hole 
as the smallest of the various radii determined above: 



-Rext — niin[i?bow, ^tidab Rmass, ^Compj ^v] • 



(13) 



From our computations, for appropriate BLR models the 
most efficient confinement mechanism turns out to be the 
bow-shock formation, thus selecting i?bow as the stellar 
envelope extension in the equation above. In this case, for 
larger distances from the star center, i.e. for R > i?bow7 
the gas physical properties are just those of the nuclear 
wind gas (i.e., n — and T — T^), since the shock 
transition is a sharp one. Using the notation of Eq. (^, the 



71* (r) 



47r rriYiR, 



■cxt'' 



(14) 



In the current framework, the stellar wind gas gener- 
ates the line emission via central radiation re-processing, 
hence the quantities i?cxt and n* determine the observed 
properties of the emission lines. An AGN wind which 
generates suitable values of the quantities n* and i?oxt 
through Eqs. (|lj) and (|l3|) can then represent the re- 
quested link between the central engine and other struc- 
tural components of the AGN itself, such as the BLR. 



3. Mass deposition in the AGN wind 

In a picture where the AGN wind expands among mass- 
losing stars, the presence of the stars will infiuence the 
AGN wind dynamics by adding mass to the wind itself. 
Two different sources of mass deposition are possible: the 
mass directly lost by the stars and that stripped by star 
collisions. This implies that nuclear wind models in such 
a framework must allow for mass deposition increasing 
the total mass flux along the radial coordinate. Indeed, 
this possibility has been taken into account in Paper I, in 
order to maintain the mass density in the wind at non- 
negligible levels in the outer region, by accounting for a 
possible mass source along the wind. In this case the total 
mass flux (mass/time, dimensionally) is not a constant, 
but is rather a function of the radial coordinate r, namely 
Mw = 47rmHnwi'w7'^ = 47rA(r) = 47r^o/m(f')- In this ex- 
pression, A{r) and the constant are dimensionally mass 
per unit time per unit solid angle, whereas f-m{f) is a di- 
mensionless function of radial distance from the central 
black hole (see Paper I, Sect. 3 for details). The constant 
Aq is deflned as Aq = ricmur'^Csc, where Csc is the sonic 
speed at the critical point Tc and ric is a characteristic 
value of the wind number density, which is related to the 
number density value at the critical point by the follow- 
ing relation: nw(''c) = n-cfmifc)- With the notation defined 
above, the mass source term [mass/time/ volume/solid an- 
gle)] is V • (mnriwVw) = {dA/dr)/r'^. 

In order to evaluate the amount of mass fed into the 
wind by stellar mass loss and by star collisions, we have 
computed the contribution from each process as described 
in the following. The total contribution, up to a distance r 
from the central black hole, to the mass flux in the nuclear 
wind from these external sources is 



Mi 



,„p = 47r[^(r) -^(ro)]. 



(15) 



where rg is the radius at the base of the wind, and, conse- 
quently, 47rA(ro) is the "intrinsic" mass flux in the nuclear 
wind, at its starting point. 

Before proceeding, we need to define a schematic de- 
scription of the spherically symmetric compact dense stel- 
lar cluster that we suppose to be located at typical BLR 
distances from the central black hole of the AGN. In the 
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present work, we assume to represent it in a strongly sim- 
plified way, that is considering just three different stellar 
components for it, main sequence stars (MSs), red giants 
(RGs), and supergiants (SGs). These three stellar compo- 
nents lose mass to the external interstellar medium, which 
is here the nuclear wind plasma, both due to their own 
stellar winds, whose relevance depends of course on the 
specific stellar component, and because of the possible 
collision-like encounters that the stars undergo. We also 
define the relative percentage of RGs and SGs, with re- 
spect to the total stellar number density of the cluster, p,, 
as /rg and /sg respectively, so that (p*)rg = /rgP* and 
(yO*)sG = fsGP*, whereas, for reasonably small values of 
the fractions /rg and /sg: we can take (p*)ms — P*^ '-e-, 
/ms — 1. In the computation of our models, we have usu- 
ally assumed a "standard" value for the relative fractions 
of RG and SG stars, that is /rg = 0.01 and /sg lO"'' 
(see ANl), but we keep those fractions generically as pa- 
rameters of the problem in what follows. 

As for the stellar mass loss from the stars of a central 
cluster characterized by a stellar number density (r) and 
a composition as defined above, the total contribution to 
the AGN nuclear wind mass loading rate, up to a radial 
distance r can be expressed as 



)ms/ms + (^4)rg/i 



RG 



(16) 

(A;f*)sG/sG| pA-^r^dr, 



where (M*)i represents the mass loss rate in the wind of 
the "i"th stellar component of the cluster. For reasons 
that will be extensively discussed in the next sections and 
are related to the very specific and extreme environment 
in which the ce ntral stellar cluster is located (see, for ex- 
ample, Taylor 1999 ), here we allow for a possible non- 
negligible contribution from main sequence (MS) stellar 
winds. 

As already mentioned, the other contribution to mass 
deposition into the AGN nuclear wind comes from stel- 
lar collisions. Following ANl and Murphy et al. (1991), 



with (i?*)i the stellar radius of stars of type "i" , and /i the 
fraction of stars of type "i" , that is assumed to be equal to 
unity for the case "i"= MS, that is for the main sequence 
star component. 

For the definition given above of Minp(r), the following 
relationship must hold 



M„p(r) = M,wind(r-) + M*coii{r). 



(19) 



Indeed, we are interested in the equation that can be ob- 
tained from the previous one by deriving it with respect to 
r and dividing by 47r, and that represents the equation for 
the sum of the mass sources from stellar matter expulsion 
processes at distance r, giving the total mass source func- 
tion per unit solid angle {dA{r) / dr) / = {Ao/r'^)df^/dr 
for the nuclear wind (see Paper I). This equation turns 
out to be the following: 



ar 

where the function F{r) is defined as follows 



(20) 



F{r) = [(M*)ms/ms + (A^*)rg/rg + (M,)sg/sgJ P*{r) + 
0.1MQ7TV{r)pl{r) x 

[(-R*)ms/ms + (^*)rg/rg + (^^»)sg/sg] • (21) 



we only take into account collisions between stars of the 
same type, thus neglecting the contribution due to MS- 
RG and MS-SG star collisions, which are characterized by 
a much smaller fractional mass loss with respect to colli- 
sions between same type stars. In the hypothesis that for 
each collision a fraction /c (that we tak e to be 0.1 follow- 
ing ANl and Begelman & Sikora 1992 ) of one solar mass 
of matter is stripped from the colliding stars and fed into 
the wind, the total mass deposition rate due to collisions 
occurred up to radial distance r, M^coiiif), is 

M*coll(r) ~ 

0.1 [(tcoii)ms(?') + (tcou)rg(»') + (Tcon)sG(?')] Mo/yr, (17) 4.1. Free parameters in the model 



From Eqs.(|20D and (|2l|), given p*{r), it is evident that for 
a given choice of the physical parameters of each stellar 
component, A{r) can be derived by solving the equation 
itself. As it will be discussed in a following Section, the 
prescription of /0*(r) is constrained by several factors, but 
p*{r) must be chosen so as to obtain a mass source func- 
tion with a sufficiently smooth behaviour. This is relevant 
for the successful integration of the nuclear wind equa- 
tions. Of course, the resulting mass source function de- 
pends on the mass loss rates that we adopt for our model 
cluster stars, as well as on the chosen values of (-R*)i, and 
of the population fractions /i, for the three stellar com- 
ponents of the cluster; this is also extensively discussed in 
the following Sections. 

A numerical comparison of the two terms on the right 
hand side of Eq. (pl|) can show whether one of the two 
processes of mass deposition in the nuclear wind predom- 
inates and, in that case, which one it is, as a function of 
r. For a given p,,{r), the result of this analysis strongly 
depends on the choice of the mass loss rates of the various 
types of stars in the model cluster. 



4. Model parameters 



where (Tcoii)i is the star collision rate (here expressed as 
number of collisions per year) for the stellar component 



defined as: 



(Tcoii)i(r) =4^'(i?,)? / {fip.)\^V{r)dr 



(18) 



As it is apparent from the discussions in Sects. 2 and 3, in 
order to derive the effective stellar envelope extension and 
its density through Eq.s(p^) and (p^, the wind model of 
Paper I must be extended to take into account the mass 
deposition function (as derived from Eqs. (|20|) and (pl|)) 
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and to include the computation of i?oxt(f ) and n*(r). This 
process implies the introduction of a number of parameters 
that can be divided in three classes, each associated with 
a different aspect of the problem. These parameters are 
summarized in the following for the sake of clarity. As 
described below, some of these parameters have almost 
no influence on the final results, some have been fixed to 
standard values and some other represent the control keys 
to variations in the outcome of our model. 

a) Parameters related to the AGN nucleus and its 
wind: central radiation field luminosity, L, characteris- 
tic number density, Uc, and temperature, T^, of the wind 
plasma at the critical point (that is at the distance at 
which the wind outflow changes its character from sub- 
sonic to supersonic). For a given luminosity the choice of 
Tc does not affect the results since a different temperature 
value essentially shifts the critical point position leaving 
all physical quantities profiles almost unaffected. We have 
shown in Paper I that the characteristic wind number den- 
sity, ric (related to the number density at the critical point, 
as defined in Sect. 3) is a critical parameter, and the sys- 
tem is very sensitive to its variations. Indeed, a solution for 
the nuclear wind equations is possible and characterized 
by a total electron scattering optical depth of the wind 
plasma ^ I (see Paper I) only over a rather narrow range 
of values for n^. Variations of within this range could 
in principle affect the BLR model and consequently might 
change the emission line properties; we discuss this issue 
in Sect. 5. 

h ) Parameters characterizing the nature of the stellar 
cluster orbiting around the nucleus: star number density 
distribution p*{r), mass-loss rate and radius (M*,i?») for 
each type of star (SG, RG, MS), and percentage of each 
type of stars present in the cluster (/j for stellar type 
"i"). These parameters enter the computation of stellar 
envelope properties, but also contribute to determine the 
mass deposition function fm{r) and hence influence the 
nuclear wind integration (see Sect. 3). Some parameters 
have been fixed to their "standard" values. This is the 
case for the star percentages and radii which are set as 
fsG = 10-^ /rg = 10-^ /ms 1 and (i?OsG = lOOi?0, 
(^*)rg = lOi?0, (i?*)MS = IRq- For others, namely /9*(r) 
and (A/*)i, different possibilities have been analysed. In 
particular for stellar mass loss rates the assumed values 
and their influence on the flnal results are described in 
the following sections. 

As far as the star number density distribution is con- 
cerned, there is no observational information about its de- 
tailed profile in regions so close to the central nucleus, 
the only possible indications come from theoretical work. 
Therefore, we have chosen /0*(r) profiles in accordance 
with two main suggestions derived from the results of 
David et al. (1987a, 1987b) analyses of a star cluster evo- 
lution around a black hole. The first suggestion regards 
the peak value for the star number density in an evolved 
cluster, that turns out to be around 10^^ stars/pc'^. The 
second one is that an evolved star distribution is concen- 
trated in a range of distances from 10^'* cm up to 10^^ cm 




16 17 
log[r(cm)] 



1 2 3 4 5 6 

log(r/rJ 

Fig. 1. Star number density as a function of the distance 
r from the AGN nucleus. We also show the abscissa axis 
for the corresponding dimensionless distance x — r/rg for 
the case of a central black hole of mass 5.6 X lO'^ Mq, 
that we have chosen to associate to a central luminosity 
L = 10^** erg/s (see Paper I). Label meaning is explained 
in the text. 



from the central black hole and steeply decreases both 
for shorter and for longer distances. We have analysed 
several different p*{r) profiles obeying these criteria, and 
present here the results for the stellar distribution profiles 
that seem to best represent our working scenario. The se- 
lected curves are shown in Fig.^ In this figure the curve 
labeled "1" is representative of the class of star density 
distribution which gives the best results for our model. 
All the results presented in the following sections refer to 
this particular choice of (r) . The other two curves have 
been illustrated since an analysis of the consequent results 
allows us to highlight separately the effects of increasing 
(r) in the external region (curve "2" ) and those of de- 
creasing p*(r) in the intermediate region (curve "3"). 

Possible differences in the stellar distribution profiles 
in the region r < 2 x 10^^ cm have proved to be unim- 
portant in the present context, owing to the fact that our 
integration of the nuclear wind equations stops at around 
that distance. We elected to retain the same profiles p*(r) 
as shown in Fig. ^, when analyzing configurations with 
different central black hole masses (that we associate with 
a different central luminosity, see Paper I), thus imply- 
ing different values for rg, and, consequently different re- 
sulting profiles as functions of the dimensionless distance 
X = r/rg. 

c) Parameters defining the physical properties of the 
stellar envelopes: stellar wind velocity profile (see Eq.(^) 
and temperature value (woi, ai,(T*)i), for each type of 
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star SG, RG, MS). This group of parameters only 

influences the properties of the line-emitting regions. 

In summary, once we have chosen the luminosity of 
the AGN for which we want to build up a model, the free 
parameters we can work with turn out to be the following 
quantities: Uc, p*{r), {M^,)i, VQi, Ui and (T,)^. The "obser- 
vational requirements" on the BLR described in the fol- 
lowing subsection dictate further limitations to the num- 
ber of parameters that significantly influence the model. 



4.2. BLR parameters and observational requirements 

Not every combination of the above listed free parame- 
ters can reproduce the observed BLR general properties. 
Observative and interpretative work has shown that def- 
inite ranges of parameters and/or quite specific features 
are common to the BLRs of most AGNs. In this section 
these observational requirements are "translated" in terms 
of our model variables. 

From the observational and interpretative points of 
view, our spherical stellar atmospheres can contribute to 
the BLR emission only if Rcxt{f) and satisfy the 

following conditions: 

a) RcKt{r) > Ai? = 10^^J7/n* cm, since photoioniza- 
tion equilibrium requires that emission lines are generated 
in shells of thickness AR, which is at least AR — 10^^C//n, 



cm (see eg. Peterson 1997). The definition of an emitting 
shell then allows an estimate of the mean density of the 
emitting plasma in the shell itself. Since the density inside 
the stellar expanding envelope varies as shown in Eq.(^, 
the mean density can be evaluated as 



Air mnvoRl 




R 



-{a+2) 



dR 



-AB. \R* J Ai?' ^ 



where both the integration limits and possibly the mass 
loss rate depend on r. In terms of the quantity just defined, 
the other conditions a BLR emitting plasma must satisfy 
read: 

b) > 10^ cm~^, since broad forbidden lines are 



not observed (Netzer 1990, Peterson 1997); 

c) if the line width is due to thermal effects only, it 
must be n*(r) < 1 x 10^^ cm~'^ to avoid line thermal- 
ization, i.e., for line emissio n to b e effectively signi ficant 
(see Rees, Netze r fc F erland |l989i Kaspi & Netzer [1999 , 
Korista & Goad |200C|) . 

In our model -Roxt('') and ^^(r) will attain different val- 
ues at different distances from the central black hole, hence 
conditions a), b) and c) will be satisfied only in a specific 
range of distances. We call ri and r2 the distances which 
delimit this interval. Stellar envelopes satisfying points a) 
to c) contribute to build up the total covering factor. Here 
we define 



air) 



TT{RcKt)lfiP*dr, 



(23) 



as the covering factor accumulated at distance r from the 
central black hole due to contributions from star type 'i" 



("i"= SG, RG, MS), obtained using Rc^t{r) computed for 
that specific star type. To obtain the total contribution to 
the covering factor due to a given stellar type, in Eq.(p3|) 
we have to integrate up to r = r2, i.e. up to the external 
boundary of the interval defining the broad line contribut- 
ing region for each stellar type, that is we have to compute 
Ci(r2) for "i"= MS,RG,SG. Finally, we define the total 
covering factor Ctot for a given BLR model as the sum of 
the total contributions, that is 

Ctot = Cms ('"2) + Crg(»'2) + Csg(»'2). 

Our model can now be compared with observations, 
testing the following points: 

1) reasonable values should be recovered for the cov- 
ering factor Ctot, w hose e stimates range from about 0.05 
up to 0.25 (Baldwin |l997| ); 

2) the typical distance characterizing the BLR, namely 
''BLR (see Sect. 1) must be contained in the intervals 
[''i,''2]i relative to the stellar types contributing to the 
definition of the value of Ctot! 

3) the covering factor relative to "broad" forbidden 
emission lines, Cforbi must be vanishing. In our model, we 
define 



Cforb = y^(Cforb)i = X! / ■^{Rcxt)ffiP*dr, 



(24) 



where 



(4 



(i?cxt)i if n4(i?ext)i] < 10« cm-3 
if n4(i?ext)i] > 10« cm-3 



(see Peterson 1997, Krolik 1999| ), therefore, a contribution 
to Cforb comes only from those envelopes whose photoion- 
ized shell interior to the confining bow shock has a density 
appropriate for forbidden line emission. 

4) the ionization parameter should be in in the range 
0.01-1. Assuming < hi' >= 2.7 Ryd as in ANl, and with 
L in erg/sec, n* in cm~^ and r in cm, the expression for 
U shown in the Sect. 1 becomes 



U{r) = 4.5 X 10" 



(25) 



The above general requirements allow us to evaluate 
the validity of our model and to restrict the range in 
which free parameters described in the above subsection 
can vary. 



5. Model results 

In Paper I, we have defined a model for an AGN nu- 
clear wind, that turned out to be essentially determined 
by the choice of the central radiation field luminosity, L, of 
ric, the characteristic number density of the wind plasma 
(see Sect. 4.1), and of the mass deposition function fm{r) 
(see Sect. 3 and Paper I). As it is discussed in Sect. 3, the 
mass source function (and therefore ultimately fmir)) for 
the nuclear wind depends on the choice of both the star 
number density distribution of the central cluster, /c*(r), 
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and the star cluster composition and properties (such as 
mass loss rates, radii and percentages of the three differ- 
ent star types that we assume as components of the stellar 
population of the cluster; see Sect. 4.1 for details). As a 
consequence, the choice of all the cluster parameters men- 
tioned above strongly influences the general properties of 
the solution for the nuclear wind as well. These, in turn, 
enter the determination of the physical quantities charac- 
terizing the stellar envelopes that we consider to be the 
line emitting sources. 

The above described complex interplay between star 
and interstellar medium properties makes it difficult to 
predict our model results. While it is obvious that the 
validity of our model could be proved only by the de- 
tailed comparison of predicted and observed line profiles, 
it seems of primary importance to understand which are 
the relevant parameters in our picture. Hence, due to the 
large number of parameters and the complexity of the 
model, a general analysis aimed at performing a first se- 
lection among all the present parameters is preferable, in 
order to get a deeper insight into the fiexibility of the 
model itself. Therefore, as a first step in the comparison 
between observed BLR features and our model results, 
we have analysed in detail the influence that each free pa- 
rameter has in the fulflUment of the observational require- 
ments presented in Sect. 4.2. The results of this analysis 
are presented hereafter. We refer to Sect. 4.1 for the spe- 
cific definition of the values chosen and fixed as "standard" 
for stellar radii and percentages in the cluster. 

Effects of changing the star density distribution . 
As we have mentioned in Sect. 4.1, point b), our system 
is very sensitive to the choice of p*. Indeed, the results 
obtained by using the three stellar distributions shown in 
Fig. |l| turn out to be significantly different. The "small" 
differences among profiles are sufficient to exclude the 
star distribution labeled with "2". In fact, any stellar dis- 
tribution (r) that is more extended in the "external" re- 
gion (i.e., for r ^ 10^* cm) than the one characterized by 
profile "1") can reproduce a covering factor in accordance 
with point 1) of Sect. 4.2, but also shows a resulting non- 
negligible forbidden line covering factor, so that it cannot 
satisfy point 3) of Sect. 4.2. On the other hand, for profiles 
similar to the one labeled with "3" , i.e. characterized by 
a smaller amount of stars, the forbidden covering factor is 
negligible, but Ctot turns out to be too low with respect to 
the range of values inferred from observations if suitably 
high values of stellar mass loss rates are not used. This 
is apparent from Eq. (p3|) , taking into account that -R^^t 
is generally proportional to the ratio M/vq (see Eq.(M)). 
Accounting for the considerations above, in the following 
we have chosen to discuss the results pertaining to the 
the star number density profile labeled with "1" in Fig. |l|, 
as representative of a possible cluster star distribution for 
broad line emitting AGNs. 

Effects of changing the nuclear wind density parameter Uc- 
As shown in Paper I (Eq.s (|l^ and (|l^)), the AGN wind 
integration system depends on the quantity /,n//m, which 
can be easily derived as 
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/m Aof^{ro) + l[r^F{r)dr 

where it is f^^ — r^F(r)/Ao from Eq.(pO[). Here the param- 
eters of the star cluster define the mass source function per 
unit solid angle, that is the function F{r), through Eq.(pT|) 
and /m(?'o) is an integration constant. Once chosen a spe- 
cific model, i.e. for fixed F{r), the relative weight of the 
terms Ao/m(ro) and /^^ r'^F{r)dr defines the amount of of 
mass flux characterizing the nuclear wind at its base with 
respect to that fed by stars. 

Hence, changing ric, which enters the definition of Aq, 
changes the function /4i//m mainly in the "internal" re- 
gion (r ro) where the term Ao/m(^o) is non-negligible 
with respect the integral term. Increasing n^, the intrinsic 
mass fiux of the nuclear wind becomes more important 
with respect to the mass flux fed by stars and the elec- 
tron scattering optical depth may become too large (as 
discussed in Paper I). For r >> ro, in the region where 
the BLR is located, the integral /^^^ r^F{r)dr is the pre- 
dominant term in the mass flux definition and a change 
in the critical density value does not have a significant ef- 
fect on the final result. A numerical analysis confirms the 
present analytical discussion, therefore we conclude that, 
once ric is chosen in the restricted range allowed by the 
nuclear wind integration (as explained in Sect. 4.1 a), its 
changes do not affect the BLR properties as derived from 
our model. 

Effects of changing the velocity profile of the stellar 
wind, i.e. the value of a, in Eq.^. 

Several integrations of our model, carried on with differ- 
ent a values, have shown that the specific velocity profile 
does not have a strong influence on the results. However, 
positive a values are generally preferable, since the corre- 
sponding models are characterized by lower values of the 
forbidden line covering factor with respect to those with 
negative a values. Perhaps, the most straightforward way 
to understand this is to go along the following line of rea- 
soning. When the stellar envelope extension is delimited 
by the bow shock (which is the case under most of the 
conditions of physical interest), the density of the enve- 
lope at its boundary, n,, is determined by the local den- 
sity of the AGN nuclear wind plasma, since oc n^, as it 
can be realized by inserting Eq.(|^) in Eq.(^. This implies 
that, going farther from the AGN central black hole, n* 
decreases, like ny^{r) does. Accounting for this condition 
in Eq.(^, it is easy to see that, for a given M, i?cxt in- 
creases with r, and it increases faster for negative values 
of a. As a consequence of this last point, for negative a 
values non-negligible contributions to the covering factor 
are expected to come from regions that are far enough (in 
the computations that we present in the following sub- 
sections, this is generally for r ^ 2 — 3 x 10'*7-g) from the 
center of the AGN, so that the stellar envelope density is 
sufficiently low to be suitable for forbidden line emission. 

Effects of changing the star envelope temperature. 
Values of the temperature in the range 1 — 3 x 10'* K, ap- 
propriate for BLR emission, have been tested, obtaining 
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no significant change in the results. This fact supports our 
simplifying assumption of a constant temperature profile 
for the overall stellar expanding envelope. 

Taking into account all the above points, we conclude 
that in our BLR model only a specific class of star num- 
ber density profiles meets the observational requirements, 
while changes in the values of the parameters n^, a and 

do not have a relevant influence on the quality of the 
reproduction of BLR properties. For this reason, in the fol- 
lowing we illustrate the general results of our exploration 
of physically reasonable ranges of the parameters whose 
changes turned out from our analysis to significantly affect 
the results as for BLR global properties, namely stellar 
mass loss rates and expansion velocity of the stellar atmo- 
spheres (see end of Sect. 4.1). To do this, we have chosen 
to discuss separately three different choices for the defini- 
tion of the mass loss rates in the envelopes of the different 
stellar types composing the cluster, and to examine the 
effects of changing the expansion velocity of the envelopes 
on a significant quantity such as the resulting covering 
factor of the proposed BLR model. Before entering the 
details of the discussion, we just want to remind that all 
the results that we show and comment in the following, 
unless otherwise specified, refer to a specific choice for 

(i) the central source luminosity, L — 10''^ erg/sec, 

(ii) the black hole mass, A/bh = 5.6 x lO^Af0 (thus 
defining rg), 

(iii) the envelope temperature, T* = 2 x 10"' K (see 
Sect. 4.2), 

(iv) the expansion velocity power law exponent, a — 
0.75 (see Sect. 4.2), 

(v) the star number density, p*(r), labeled "1" in Fig. |l|, 

where choices (iii), (iv) and (v) are justified by the above 
discussion. 

As for the stellar mass loss rates, we have chosen to 
discuss 

a) the simple case of "standard" mass loss rates for 
stellar winds of stars of a given type, as known from galac- 
tic studies (see below for the details and values. Sect. 5.1); 

b) a case in which we try to account for the possibility 
of "enhanced" mass loss of MS, RG, and SG stars, due to 
the very specific context in which they are set, namely a 
strong X-ray illumination from the central source in the 
AGN (see Sect. 5.2); this would lead to an expected en- 
hancement of mass loss rates increasing with decreasing 
r, with a behaviour following the radiation flux increase 
going inwards, that is cx r~^; 

c) a third case in which we suppose the mass loss rates 
do again depend on the distance from the central black 
hole r, but with a flatter increase going inwards (Sect. 5.3). 
The physical mechanism responsible for this configuration 
is not specified in this case, but we refer to Sect. 6 for a 
brief discussion of possible processes. 



5.1. Standard mass loss rates 

The first reasonable step is that of assuming "stan- 
dard" mass loss rates for the stars in the cluster (see, for 
example, Lamers & Cassinelli 1999| ), that is 



(M, 



*)SG 



IQ-^MQ/yr, 



(M,)rg - lO-»M0/yr, 
(Af,)MS - lO-i^Mo/yr. 



ms 

With this choice for the mass loss rates of the stars in 
the central cluster, from Eq. ( pO[ ) it turns out that, in the 
inner region of the cluster (i.e., for r ^ lO^'rg), the mass 
flux input in the nuclear wind is strongly dominated by 
coUisional effects (i.e., mass loss from collisions between 
stars) . On the contrary, farther from the AGN center there 
is a vast region in which predominance of the stellar wind 
mass loss sets on. We then derive the corresponding AGN 
nuclear wind model, and, consistently, the main physical 
parameters of the BLR emitting envelopes. 

For each stellar type these are: i) the evaluation of the 
characteristic extension for the expanding envelope, i?ext 
(see Eq. (^3|)), and the physical mechanism defining it; 
ii) the average density of the envelope over the emitting 
region of the envelope itself, n,, (see Eq. (p2|)); iii) the 
resulting ionization parameter U (Eq.(p5|)); iv) the cover- 
ing factor, as a function of the distance r from the black 
hole in the interval [ri, defining the BLR (see Eq. (p3|)). 
A typical behaviour (which, we want to stress, is gener- 
ally common to the three cases for stellar mass loss rates 
that we present in this work) corresponds to an increas- 
ing envelope extension, and to a decreasing density with 
r increasing. Of course, the covering factor Ci(r) is an 
increasing function of r over the model BLR range of dis- 
tances, since the contributions from the stars in each shell 
at a given distance from the black hole accumulate in the 
integral defined by Eq. (|2^). 

First, we want to focus our attention on the results we 
obtain from the model in terms of global parameters of 
the BLR, such as the integrated covering factor, defined 
as Ci(r2) for each stellar type "i". We briefly describe the 
r dependence of the first three physical quantities men- 
tioned above later on in the present Section and in Sect. 6. 

The general outcome for this case is exemplified in 
Fig. H, showing the resulting covering factor contribu- 
tions (computed from Eq. ( p3| ) with r — for the mod- 
eled BLR, as functions of the terminal velocity Uterm = 
(i?cxt), characterizing the stellar envelopes that are as- 
sumed to contribute to the broad line emission. We also 
plot the sum of the contributions due to the different stel- 
lar types composing our model central cluster as Ctot, 
again as a function of Wtcrm- We obtain the variation of 
"yterm for cach stellar type, by varying (wo)ii which, in turn, 
also induces variations in the extension of the envelope, at 
least when the confinement mechanism for the envelope 
itself is the bow shock (see Sect. 2, Eqs. (|^) and (p^)); 
changing this parameter is therefore significant to the de- 
termination of the covering factor. To give an idea of the 
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range of {vo)i, we report here the values corresponding to 
the terminal velocities utcrm = 0.1 and 2.1 km/s, which 
are close to the extremes of the terminal velocity range 
plotted in Fig. 2 and representative of "low" and "high" 
terminal velocities respectively. In the present case, for the 
star types that do give contribution to the covering fac- 
tor, the values of vq are (wo)RG('yterm = 0.1 km/s) = 1.6 x 
10-3 km/s, (wo)sG(wtcrm = 0.1 km/s) ^ 3.7 X 10-3 km/s, 
and (uo)RG(wtcrm = 2.1 km/s) = 1.8 x 10-^ km/s, 
(i^o)sG(wterm = 2.1 km/s) = 4.2 x 10"^ km/s. 

It is important to note that the plotted value of Ctot 
is the one that we obtain when we suppose that the three 
different star type envelopes are endowed with the same 
value of the terminal velocity, which is not necessarily the 
case. In fact, one can obtain different values for the to- 
tal covering factor for the BLR model, with respect to 
those that we present in the plot, by assuming different 
^^term valucs for cach single species of stars, that is pick- 
ing contributions from the various star types at different 
wtcrm along the curves defining the specific covering factor 
for a given stellar type. However, we want to stress that, 
due to the fact that the profiles Ci(wtcrm) for the three 
stellar types show a similar behaviour, which is generally 
decreasing for increasing Vterm (apart from the very low 
velocity end of the range shown in the figures), it turns 
out that for any chosen value of Vta-m the total covering 
factor we plot in Fig. ^ is the maximum value that can be 
attained if we suppose that the different contributing stel- 
lar type envelopes are characterized by terminal velocity 
values that are possibly different, but in any case larger 
than the one chosen. On the contrary, larger Ctot values 
can be obtained by allowing for different terminal veloci- 
ties of different stellar type envelopes, provided the Vta-m 
value of at least one of the contributing types is lower than 
the one chosen for the exemplification. These are general 
considerations and they refer both to the present and to 
the other two cases we discuss in the following subsections. 

An inspection of the solutions for the BLR physical 
parameters shows that the bow shock mechanism is in 
general the most efficient in confining the stellar envelopes 
and, consequently, determines their extension for all the 
three different stellar types. As a consequence, for this 
specific choice of the stellar mass loss rates, the envelope 
radii for RG stars and SG stars turn out to scale so as to 
give essentially the same contribution to the total covering 
factor for the BLR. As it is apparent from Fig. ^, this 
is no longer true for low velocities ( ^ 0.12 km/s), when 
the contribution to the covering factor due to SG stars 
"decouples" from that of RGs. The physical reason for 
this is in the fact that at these low velocities the most 
efficient confining mechanism for SGs all through the BLR 
extension becomes tidal disruption of the envelope beyond 
a given radius (see Eq. (||)), whereas for RGs the bow 
shock mechanism is still dominant. 

To illustrate the behaviour of the BLR physical quan- 
tities and the characteristic extension of our BLR model 
(i.e., the [ri,r2] interval over which conditions a) to c) of 
Sect. 4.2 for the emitting gas structures are fulfilled), we 
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Fig. 2. The total value of the covering factor, Ctot, de- 
fined as the sum of the contributions due to the envelopes 
of different stellar types, as well as the distinct contri- 
butions to this value from the various stellar types, are 
shown as functions of the terminal velocity of expansion 
^torm (in km/s), characterizing the contributing envelopes. 
The present figure refers to the case for "standard" mass 
loss rates. Notice that the contribution due to main se- 
quence stars in this "standard" case does not even appear 
in the plot, since it is absolutely negligible with respect to 
those due to RG and SG stars. 

have chosen to plot in the two panels of Fig. |^ (each one 
referring to one of the two stellar types that give substan- 
tial contribution to the broad line formation, namely RGs 
and SGs in this case) the three quantities i?cxt, n* and U. 
These are functions of the distance from the central black 
hole r, and we plot two curves for each quantity, corre- 
sponding to the two representative values of the expan- 
sion velocity parameter Wtcim previously defined, namely 
I'tcrm = 0.1 km/s and Wtorm = 2.1 km/s, close to the ex- 
tremes of the velocity range plotted in Fig. |2[ For the sake 
of clarity, we note that the utcrm values chosen refer to 
those attained for envelopes at r = r2, that is at the ex- 
ternal border of the broad line emitting region. For each 
BLR physical parameter the curves obtained for interme- 
diate terminal velocity values typically lie in the region 
between the 0.1 and 2.1 curves. Therefore, each of the two 
panels of Fig. ^ essentially shows the ranges of values for 
the BLR parameters for the models corresponding to the 
present choice of mass loss rates for the stars of the cen- 
tral cluster. It is interesting to notice that for RGs the 
extension of the BLR contributing region is independent 
of the terminal velocity value, whereas this is not the case 
for SGs. This is again due to the fact that for RGs the 
envelope radius is defined by the bow-shock mechanism, 
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Fig. 3. Physical parameters characterizing stellar en- 
velopes contributing to the BLR as functions of the nor- 
malized distance r/rg from the central black hole; these 
are the mean density of the envelope, the estimate 

of the envelope extension, i?oxt ('') , and the ionization pa- 
rameter, U{r). The interval over which the physical quan- 
tities are plotted represents the "BLR" interval [?'i,r2] as 
defined in Sect. 4.2; the labels 0.1 and 2.1 identify the spe- 
cific value (in km/s) of the terminal velocity parameter for 
which they are computed. The upper panel refers to RGs' 
envelopes, while the lower one shows the results pertain- 
ing to the envelopes of SG stars, in the case of "standard" 
mass loss rates. 



whereas at low wtcrm it becomes (i?oxt)sG = (-Rtidai)sG- 
Indeed, it is easy to see that, when the bow-shock is the 
confining physical mechanism for the stellar envelope, and 
the expansion velocity is subsonic, n»(r) (see Eqs. (6) and 
(7)) is essentially independent of the stellar wind veloc- 
ity itself. Moreover, when the envelope gas temperature is 
chosen to be the same for the three stellar types, n*(r) is 
independent of the stellar type as well. The same consid- 
erations do hold for the average density in the emitting 
portion of the envelope, (see Eq. (|2^)), which is the 
quantity that we plot in Fig. |^, since in general it turns 
out that = n*. In fact, our computations have shown 
that the emitting shell of an envelope is characterized by 
a size AR << i?cxt, so that, performing an analytical 
integration of Eq. ( p2[ ) and using a first order approxima- 
tion in Ai?/i?oxt(<< 1)) one ends up with n» = As 
a result, the limiting conditions on for line emission 
(points b) and c) of Sect. 4.2) are met at the same dis- 
tance from the black hole, r, independently of the chosen 
value for the velocity parameter, since (see Eqs. (^ and 
(^) they only depend on the nuclear wind physical prop- 
erties, namely nw(?'), the nuclear wind density, Tw(r), its 
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Fig. 4. Temperature (dashed line) and number density 
(continuous line) profiles for the AGN nuclear wind so- 
lution obtained for the case of mass deposition defined 
by the choice of "standard" mass loss rates for the clus- 
ter stars. For comparison, we also plot as the dotted curve 
the "most extreme" density profile that we have obtained, 
namely the one resulting from the solution for the nuclear 
wind with the mass deposition defined by the choice of 
mass loss rates described as case c) in Sect. 5 and dis- 
cussed in detail in forthcoming Sect. 5.3. 



temperature profile, and on the velocity defining the ram 
pressure at the envelope boundary (see Sect. 2). We want 
to stress here that these same considerations do hold inde- 
pendently of the choice of mass loss rates, therefore they 
similarly apply to the cases discussed in the next sub- 
sections as well. For completeness, in Fig. 4 we show the 
behaviour of the AGN wind physical quantities n^{r) and 
T-v,{r), as obtained from the nuclear wind system resolu- 
tion for this specific choice of the mass deposition. For this 
density profile the resulting total optical depth to electron 
scattering for the nuclear wind plasma is tt — 0.12. Note 
that, for comparison, we have also plotted in Fig. 4 the riw 
density profile (dotted curve) resulting from the compu- 
tation of the solution for the nuclear wind with a different 
mass deposition, that is the one that will be extensively 
discussed in Sect. 5.3, and corresponds to our "most ex- 
treme" (i.e., strongest Mjnp) choice for mass deposition in 
the present work; the temperature profile corresponding 
to this latter case is not shown, since the differences with 
the dashed curve of Fig. 4 are qualitatively irrelevant. 

Going back to the limiting conditions (points b) and c) 
of Sect. 4.2) for the BLR contributing envelope density n,, 
for the present case they define the interval [ri,r2]RG — 
[237rg, 1.7 X 10\] = [3.95 x lO^^ cm,2.8 x 10^^ cm], thus 
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including any reasonable estimate of the characteristic size 
''BLR for our chosen value of L. 

From Fig. it is apparent that, within the "standard" 
mass loss rate case defined above, significant covering fac- 
tors can be recovered only when the terminal velocities 
of expansion for the stellar envelopes are quite low, as 
compared with the typical sonic speed for the envelope 
gas (for ~ 2. X 10"* K, we get Cs ~ 21 km/s). To 
obtain Ctot — 0.1, wtorm should be around 0.027 km/s. 
Indeed, this is a very low value for the expansion velocity 
of a "normal" stellar expanding atmosphere. However, we 
would like to mention that similarly low values have been 



presented in AN2 (1997) in their own attempt to model 



the BLR structure as an ensemble of expanding stellar en- 
velopes. In fact, transposing their notation to ours, these 
authors (AN2) assume a = —0.5, a decreasing function of 
distance from the star surface for the stellar wind velocity, 
and they choose to analyze models with vq ~ 0.1 km/s, 
where the value of is, in their case, the peak value of 
the stellar wind velocity. 

We want to note that the present case of "standard" 
mass loss rates seems to be not very much relevant to 
the construction of a physical BLR model for another rea- 
son, namely the fact that the calculation of the resulting 
covering factor for forbidden lines, which are actually not 
seen in the BLR, gives uncomfortable results. In fact, the 
global value of Cforb (as defined in Eq. (|4|)) turns out 
to reach a significant fraction of the obtained Ctot- For 
example, when Ctot — 0.1, we get for the present case 
Cforb — 0.037, which is not negligible. This is related to the 
characteristics of the nuclear wind solution for the present 
case, in which not much mass is input in the nuclear wind 
as a consequence of stellar mass loss, and, therefore, the 
nuclear wind plasma density is rather low in the "ex- 
ternal" (i.e., at r ^ 10 tblr) region; this allows for the 
formation of larger expanding envelopes, characterized by 
density values that are sufficiently low to contribute to the 
increase of the covering factor for forbidden lines (see the 
trends shown in Fig. ^ for the envelope physical parame- 
ters). 

As a final remark, we note that the total number of 
stars contributing to the BLR, that is essentially the num- 
ber of RGs in the interval of distances [ri, for this case, 
turns out to be 2 x lO'^. The total number of stars present 
within the bounds of the BLR for this case is therefore 
~ 2 X 10^. 

5.2. Enhanced mass loss rates 

It is by now well known that in a close binary system, 
intense X-ray illumination from the primary star can have 
noticeable effects on the secondary star. The best studied 
example is that of HZ Her (in the system HZ Her/Her X- 
1), whose illuminated face originates an "X-ray-excited" 
wind (see e.g. Hameury et al. 1993). The effect is related 
to the intensity of the X-ray flux at the stellar surface, and 
from this intensity it is possible to estimate, through an 
efficiency factor, 77, the fraction of luminosity that gets 
transformed in kinetic power of the wind. By physical 



analogy, we suppose that the stars of the AGN central 
cluster may suffer from a similar effect, since they are 
strongly illuminated by the hard radiation of the central 
radiation field of the AGN itself (see, for example, Taylor 
1999 ). We can therefore try an analogous (see Hameury 
et al. 1993) estimate of the transfer of kinetic power into 
the stellar wind of stars belonging to our central cluster 
in a radio-quiet AGN, from the X-ray luminosity, Lx, of 
the radiation field of the AGN central source, as follows: 

vLx(^) c^vlM,. (27) 



The value of 77 computed for Her X-1 is 3 x 10~ ^ and it 
must be in any case < 0.1 (Hameury et al. |1993| ). Taking 
into account that Lx < L , we assume ijLx — 10~'^L. We 
notice here that, in order to produce BLR models that 
meet the general observational requirements discussed in 
Sect. 4.2, the resulting velocities for the expanding stellar 
envelopes are significantly lower than the sonic speed (by 
an order of magnitude or even more). Therefore, it seems 
appropriate to derive an estimate of the "enhanced" mass 
loss rate for any given type of stars in the cluster with 
a value close to the ones that correspond to models 
that better fulfill the observational conditions on the BLR. 
In fact, wc choose to use ~ 3 km/s to infer an order 
of magnitude value for the induced mass loss rate of a 
star in the central cluster of the AGN from relation (p7[). 
For L — 10'*'* erg/sec, we thus find an estimate for the 
enhanced mass loss rates 



(M*) 



MS 



(Af, 



8 X 10-2(r/rg)-2MQ/yr, 



*)RG 



8(r/rg)-2MQ/yr, 



(M. 



*)SG 



X 102(r/rg)-2Mo/yr. 



(28) 



(29) 



(30) 



We have computed models using for the three types 
of stars considered the mass- loss rates given by Eqs. (28) 
to (^0|) for r such that those relations give values larger 
than the "standard" ones defined in Sect. 5.1, and main- 
taining, on the contrary, those same "standard" values 
when the equations above would determine smaller values. 
At r ^ J'BLR the above estimates would give {M*)ms — 
2.4 X lO-^Me/yr, (M*)rg - 2.4 x IQ-^Mq/jt, and 
(Af*)sG — 2.4 X 10~'*Mo/yr. These evaluations point out 
an interesting consequence of this enhancement mecha- 
nism for stellar mass loss rates. In fact, at distances com- 
parable to the characteristic BLR radius, this process does 
change substantially the "intrinsic" red giant and super- 
giant mass loss rates, by a factor ^ 10^, but, in the spe- 
cific environment considered here, it has much more influ- 
ence on the effective mass loss rate of main sequence stars, 
that turns out to be increased by more than a factor 10^ 
at these distances, relevant to the BLR. Therefore, the 
first important difference with respect to the "standard" 
mass loss rate case presented in Sect. 5.1 is that the con- 
tribution of MS stars to the total covering factor is no 
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Fig. 5. The total value of the covering factor, Ctot, de- 
fined as the sum of the contributions due to the envelopes 
of different stellar types, as well as the distinct contri- 
butions to this value from the various stellar types, are 
shown as functions of the terminal velocity of expansion 
^^tcrm (in km/s), characterizing the contributing envelopes. 
The present figure refers to the case for "enhanced" mass 
loss rates, and, differently from the case shown in Fig. ||, 
the contribution due to main sequence stars is quite signif- 
icant, and actually predominant with respect to RG and 
SG contributions. 

longer negligible for any value of the expansion velocity of 
the emitting envelopes. This is immediately apparent at a 
glance from Fig. ^, which is the analogous to Fig. ^ dis- 
cussed in the previous Sect. 5.1 for the case of "standard" 
mass loss rates. The plot in Fig. ^ shows that the strong 
increase in mass loss rate for MS stars produces as a con- 
sequence a significant contribution to the covering factor. 
In fact, the extension of the emitting stellar envelopes of 
MS stars is determined by the bow shock mechanism for 
basically any expansion velocity and the envelope radius 
is determined by Eq. (j?]); therefore, main sequence star 
envelopes turn out to be much larger than in the case of 
Sect. 5.1, due to the strong increase in {M^)ms- Again, the 
number density n^{r) and temperature TW(r) of the nu- 
clear wind are computed as the consistent solution of the 
nuclear wind system of equations with the present mass 
deposition choice. In this case, the total optical depth to 
electron scattering for the nuclear wind plasma turns out 
to be tt — 0.15. 

It is interesting to note that, due to the decrease in 
{M^)i moving outwards from the central radiation source, 
the MS star envelope extension, Rcxt, shows an extremely 
slow increase (if not a similarly extremely slow decrease 
locally) with increasing r. This can be seen directly from 
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Fig. 6. Physical parameters characterizing stellar en- 
velopes contributing to the BLR as functions of the nor- 
malized distance r/rg from the central black hole; these 
are the mean density of the envelope, n*{r), the estimate 
of the envelope extension, i?ext(?'), and the ionization pa- 
rameter, U{r). The interval over which the physical quan- 
tities are plotted represents the "BLR" interval [ri,r2] as 
defined in Sect. 4.2; the labels 0.1 and 2.1 identify the spe- 
cific value (in km/s) of the terminal velocity parameter for 
which they are computed. The upper panel refers to main 
sequence star envelopes, while the lower one shows the re- 
sults pertaining to the envelopes of RG stars, in the case 
of "enhanced" mass loss rates. 

Fig. ^, similar to Fig. ^ for the "standard" mass loss case 
discussed in Sect. 5.1 (and we refer to that Section for a 
detailed explanation of the notation and of the choice of 
the curves plotted); only, in the present case, we show the 
results referring to MS stars and RG stars in the upper and 
lower panel respectively, since these two types of stars are 
the ones that mostly contribute to the present BLR model. 
As for RG envelopes, their extension is mostly bow shock 
determined, apart from the very inner shells of the BLR, in 
which the envelope radius is "tidally" defined; this "tidally 
dominated" portion of the RGs contribution to the BLR 
slowly extends to larger r distances with decreasing the 
expansion velocity parameter. We note that this does not 
influence the behaviour of the curve CRG(i'tcrm), since the 
portion of BLR over which the envelope radius turns out 
to be "tidal" is in any case not much extended in r and, 
especially, it is restricted to the inner region of the BLR 
contributing interval, whereas most of the contribution 
to the integrated covering factor Crg(''2) accumulates at 
larger values of the distance r, in the "bow shock domi- 
nated" portion of the interval [?'i,?'2]rg- Analogously to 
the MS envelopes' case, in the "bow shock dominated" 
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region the envelope radius of RG expanding atmospheres 
is more or less constant with r, and this can be seen from 
the lower panel of Fig. ^. The resulting interval of dis- 
tances defining the extension of the MS stars contribu- 
tion to the BLR is [ri,r2]MS =^ [239rg,2.2 x lO^'rg] = 
[4. X 10^^cm,3.7 x lO^^cm]. RG stars contribute in a region 
[''ij ''2]rg whose extension shows a dependence on the ve- 
locity parameter Utcrm, as it can be seen from Fig. ||, but 
it does anyway contain the estimated value of tblr, inde- 
pendently of titerm itself. As for the mean envelope den- 
sity ri*, for the conditions in which the envelope radius is 
-Rcxt — -Rbow, we refer to the considerations presented in 
Sect. 5.1, that still hold. 

From Fig. || it is apparent that, at any value of the 
expansion velocity, the contribution to the total cover- 
ing factor due to MS stars is dominant, but comparable 
in order of magnitude with the RG contribution. On the 
contrary, the contribution of SG envelopes turns out to be 
much less significant in this case. An analysis of the SG 
envelope extension as a function of r shows that when the 
velocity parameter wtcrm becomes very small (< 0.2 km/s), 
the dominant confinement mechanism for the envelopes 
becomes generally different from the bow-shock one, sim- 
ilarly to what was discussed for the "standard" mass loss 
case in Sect. 5.1, thus leading to the change in slope of 
the curve Csg — C'sG(wtorm) observed in Fig. ^. As for the 
definition of the range of values for the velocity param- 
eter (uo)i for the present case, we report representative 
values that are analogous to those defined in Sect. 5.1 for 
the "standard" mass loss rate case, namely (wo)Ms('f^tcrm — 
0.1 km/s) = 1.4 x 10"3 km/s, (t)o)RG(wtorm = 0.1 km/s) = 
1.4 x 10^3 km/s, (t;o)sG(i'torm = 0.1 km/s) = 1.4 x 
10"^ km/s, and (^^o)Ms(^^term = 2.1 km/s) = 8.8 x 
10-2 km/s, (fo)RG(t'term = 2.1 km/s) = 8.8 X 10-2 km/s, 
{vQ)sG{vtcmi — 2.1 km/s) — 8.8 x 10^^ km/s. Notice 
that here we have given the values for main sequence 
stars as well, since they do contribute to the covering 
factor in this specific case. Also, the values we report 
are exactly the same independently of the stellar type 
and this is due to the specific choice of the mass loss 
rates for the three stellar types, which turn out to give 
(M^/RDms = {M^/RI)kg = {M-,/Rl)sG, thus resulting 
in the same {Rcxt/R*) value at the BLR external border, 
defined by the condition n» = 10* cm^^^ for any stellar 
type (see Eq. (|)). 

It is easy to see that in the present case, we can build 
models corresponding to values around 0.1 for the to- 
tal covering factor for envelope expansion velocities sig- 
nificantly larger than what results in the case of "stan- 
dard" mass loss rates. In fact, we obtain Ctot — 0.1 for 
I'term ~ 0.4 km/s, which is more than one order of magni- 
tude larger than the result obtained in the case discussed 
in Sect. 5.1. However, this velocities still are in a strongly 
subsonic range of values. Another improvement, with re- 
spect to the global results of the first case (discussed in 
Sect. 5.1), is the fact that for the present choice of stellar 
mass loss rates, the ratio Cforb/C'tot turns out to be quite 
smaller, typically ranging between 0.02 and 0.04 for the 



range of velocity parameter shown in Fig. || thus setting 
the value of Cforb significantly close to "vanishing" . 

Finally, we want to stress another important differ- 
ence with respect to the case of the standard mass loss 
model, namely the fact that in the present scenario we 
have a much larger number of distinct emitting units con- 
stituting the BLR structure, since we can include all main 
sequence stars of the central cluster characterized by dis- 
tances r from the central black hole that fall in the BLR 
contributing interval ['ri,r2]MS- For example, for models 
giving Ctot — 0.1, we get an estimate of the total emitting 
units in the BLR (iV,)tot ^ 2.7 x lO"^. 



5.3. "Intermediate" mass loss rates 

As a third case for stellar mass loss rates, we have cho- 
sen another set of r-dependent (M*);. In this case, how- 
ever, we do not specify the physical motivation lying be- 
neath the choice, like in Sect. 5.2, and we only rely on the 
fact that other mechanisms could be responsible of gen- 
erating a condition in which mass loss rates of stars in 
the central cluster of an AGN could be altered (and en- 
hanced) with respect to the "standard" values (see Sect. 6 
for a brief discussion). This possibility has been taken into 
account also in previous analyses of BLR models based 
on stellar envelopes as individual contributors to the line 
emission. We refer, for instance, to the work of Scoville 



and Norman (1988), in which the authors consider giant 
star envelopes in a central cluster as emitting units; they 
suppose that for the giant stars an "induced" mass loss, 
that turns out to depend on the distance r, is generated by 
some unspecified external heating from the central source. 
Therefore, our choice in the present Section is similar to 



the one of Scoville and Norman (1988), since we do not 



analyse the physical mechanism responsible for affecting 
the stellar mass loss rates. 

We have thus parameterized the mass loss rates that 
we take into account in the present Section as follows: 

10 f T 



(M* 



MS 



9 X 10" 



(M*)rg = 3 X 10- 



V 10*^7 



-1 ( ^ 
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V 10^7 
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Mo/yr, 



(M,)sG - 3 X 10- 



10^7 



-0 



(31) 

(32) 
(33) 



where we have chosen [3 = 0.3. For any type of star, the 
equations above define the mass loss rate that we use for 
the models we discuss in this Section for r up to the dis- 
tance at which the value of the mass loss rate itself (which 
is a decreasing function of r) reaches the "standard" value 
mentioned at the beginning of Sect. 5.1, whereas for larger 
distances we discard Eqs. (^)-(|33|) and we assume that 
the mass loss rate is constant and given by the "standard" 
values of Sect. 5.1. As for the nuclear wind number density 
and temperature profiles (nw(7') and Tvj{r)) derived from 
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the nuclear wind solution referring to the present case for 
mass deposition, we remind that n^rir) is shown in Fig. 4 
as the dotted curve. The nuclear wind temperature profile 
shows no significant differences with respect to the one 
shown in Fig. 4, from a qualitative point of view. The re- 
sulting total optical depth to electron scattering for the 
nuclear wind plasma is here tt — 0.4. 

We want to note that, with respect to the case pre- 
sented in the previous Section, for red giants the present 
choice of mass loss rate function has significantly larger 
values over most of the region from which we expect to 
have BLR contributions. A similar condition holds for MS 
stars as well, but the effect is weaker than for RGs. As for 
SG stars, the situation is less well defined, but (Af*)sG(^) 
is effectively larger than what we get in Sect. 5.2 over 
the outer portion of the region contributing to the BLR 
emission for SGs. As a consequence of this, the resulting 
contributions of the different types of stars to the total 
covering factor, that we again show as functions of the 
expansion velocity parameter Wtcrm in Fig- 0; turn out to 
have a different relative importance with respect to the 
"enhanced" mass loss rate case of Sect. 5.2. In fact, in the 
present case RG stars are the dominant contributors to 
the total covering factor, but the contribution of main se- 
quence stars is significant (~ 22% of Ctot) at any value of 
the velocity parameter Wtcrm- On the contrary, the contri- 
bution due to much less numerous SG stellar component 
does not increase as much, and we get a Csg estimate 
which is < 7% of the total value of the covering factor for 
any ?;tcrm- 

A closer inspection of the BLR models we have used 
to build the curves in Fig. |^ shows again that the bow- 
shock mechanism for confinement of the stellar expanding 
envelopes is the one that effectively determines the exten- 
sion of the envelopes under any condition for MS stars. 
For RG envelopes this is still true, except for the very low 
velocity end of the interval shown, in which case other 
confinement mechanisms dominate in the inner portion 
(r £ 500rg) of the BLR contributing region, [ri, r2]RG- On 
the other hand, SG envelopes are confined basically by 
bow-shock mechanism for Vterm > 0.3 km/s, whereas for 
lower values of the velocity parameter other confinement 
mechanisms become more efficient in defining the envelope 
extension, just like the case analysed in Sect. 5.2. 

Regarding the definition of the range of values for 
the velocity parameter (uo)i for the present case, again 
we report representative values analogous to those de- 
fined in Sects. 5.1 and 5.2 for the previous mass 
loss rate cases, namely (wo)Ms(wtcrm = 0.1 km/s) = 
5.8 X IQ--^ km/s, (wo)RG(wtorm = 0.1 km/s) = 3.7 x 
10-"* km/s, (wo)sG(wtorm = 0.1 km/s) = 8.8 x 
10^'' km/s, and (wo)Ms(^^tcrm — 2.1 km/s) = 3.4 x 
10-2 km/s, (?;o)RG(?^torm = 2.1 km/s) = 2.1 x 10-2 km/s, 
('^o)sG(i^torm = 2.1 km/s) — 5.1 X 10-2 km/s. We have 
given the values for main sequence stars as well, since they 
give a substantial contribution to the covering factor also 
in this case. 
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Fig. 7. The total value of the covering factor, Ctot, de- 
fined as the sum of the contributions due to the envelopes 
of different stellar types, as well as the distinct contribu- 
tions to this value from the various stellar types, are shown 
as functions of the terminal velocity of expansion Utcrm (in 
km/s), characterizing the contributing envelopes. In the 
present "intermediate" case, both RGs and MSs give sub- 
stantial contribution to the definition of the total covering 
factor 

Similarly to Figs. ^ and [g] for the previous two cases. 
Fig. H shows the profiles of the relevant physical parame- 
ters for stellar envelopes contributing to the BLR as func- 
tions of the distance r. We show the results pertaining 
to the same two representative values of the velocity pa- 
rameter vtcrm, namely 0.1 km/s and 2.1 km/s. The exten- 
sion of the BLR contributing interval for MS envelopes 
is [ri,r2]MS = [410rg,3.2 x lO^rg] ~ [6.8 x lO^^ cm, 
5.3 X 10^^ cm], independently of the velocity parameter, 
again because of the fact that the extension of the en- 
velopes is always defined by the bow-shock mechanism. 
For RG envelopes, the BLR contributing interval has the 
same upper limit r2 as that for MS contributors, whereas 
the inner boundary ri shows a dependence on the velocity 
parameter (see previous paragraph). 

As it is easy to see from Fig. |^, the present choice of 
mass loss rates for the central cluster stars turns out to 
be the most favourable one from the point of view of the 
value of the covering factor that can be recovered with 
our models. In particular, the interval of observationally 
inferred values for a Seyfert 1 BLR, defined in Sect. 4.2 to 
be ranging from 0.05 up to around 0.25, can be obtained 
with stellar envelopes characterized by expansion veloc- 
ities in the range 0.5 — 3.0 km/s, which corresponds to 
subsonic, but physically sensible values for an expanding 
stellar atmosphere. As for the values of Cforb, they turn 
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Fig. 8. Physical parameters characterizing stellar en- 
velopes contributing to the BLR as functions of the nor- 
malized distance r/rg from the central black hole; these 
are the mean density of the envelope, the estimate 

of the envelope extension, i?oxt ('') , and the ionization pa- 
rameter, U{r). The interval over which the physical quan- 
tities are plotted represents the "BLR" interval [?'i,r2] as 
defined in Sect. 4.2; the labels 0.1 and 2.1 identify the spe- 
cific value (in km/s) of the terminal velocity parameter for 
which they are computed. We show the results for both 
MSs (upper panel) and RGs (lower panel) , since they both 
significantly contribute to the total covering factor in the 
present case of "intermediate" mass loss rates. 

out to be sufficiently small to be considered "vanishing" , 
similar to what happens in the case for "enhanced" mass 
loss rates of Sect. 5.2, since they at most reach 6% of Ctot- 
As a final consideration for the intermediate mass 
loss scenario we estimate total number of emitting units. 
Indeed, similar to the case for "enhanced" mass loss rates, 
we can include MS star envelopes in the total number of 
individual contributors to the model BLR, since their con- 
tribution to the total covering factor is also substantial in 
this case. This allows an estimate of (iV*)tot ~ 3.2 x 10^ 
stars; this number is even a little higher that the one es- 
timated for the case of "enhanced" mass loss rates; this is 
due to the fact that in the present case the external bound- 
ary of the BLR contributing interval of distances is larger 
than that of the case discussed in Sect. 5.2, thus allowing 
to include a larger number of stars in the computation. 



6. General Discussion 

The results presented in the previous section show that 
our model is capable of reproducing the observational re- 



quirements listed in Sect. 4.2 (points 1-4). In fact, from 
Fig.s H, H and ^ it is evident that the ionization parameter 
U is in the expected range of values, while tbrl is con- 
tained in the intervals [ri,r2] where the conditions a)-d) 
of Sect. 4.2 are satisfied. As far as the covering factor is 
concerned, it is evident that its value strongly depends on 
the choice adopted for the stellar mass loss rates (compare 
Fig. |, I and 0). In the case of enhanced and intermediate 
mass loss rates, values in the range 0.05- 0.25 can be re- 
covered by choosing suitable values of the stellar envelope 
expansion velocities, while for standard mass loss rates 
this is possible only for very low velocity values. Hence, 
our values of Ctot tend to exclude stellar envelopes with 
standard mass loss rates as plausible emitting units for a 
BLR model. On the other hand, regarding the covering 
factor for "broad" forbidden line emission, we recall that 
the star density distribution we have adopted has been 
chosen so that in all the presented cases Cforb is mini- 
mized (see Sect. 5); in fact, it turns out to be neghgible 
in both the cases presented in Sects. 5.2 and 5.3. We re- 
call here that in the "standard" mass loss rate case, Ctorb 
is not negligible and we consider this one of the reasons 
why that particular model is not viable for building up a 
realistic BLR model. 

While the analysis of this paper has produced strong 
constraints to most of the free model parameters, for some, 
namely the stellar distribution, the mass loss rates and 
the expansion velocity, the model cannot operate a defi- 
nite choice. Any further test of the model and selection of 
these parameters requires more information coming from 
the comparison between the observed line profiles and rel- 
ative intensities, and those predicted by the model. As 
we have already mentioned, we have planned to compute 
the detailed line profiles in a forthcoming paper, having 
now obtained a deeper insight in the role that different 
parameters play in the model. However, at this stage we 
would like to compare our model with th e results of some 
other authors (Goad fc Ko ratkar 1998 , Kaspi & Netzer 
Goad I 



1999| , Korista & Goad gOOOj) who made attempts to infer 
from line profiles and intensities some "average" physical 
conditions in BLR condensations. 

In the last years much work has been devoted to in- 
fer the BLR physical configuration starting from the ob- 
served BLR properties. Since the BLR emitting gas has 



been shown to be stratified in density (see Peterson 1993 



and Netzer & Peterson 1997 for reviews) the idea was 
to derive both the dependence of the characteristic num- 
ber density of the gas in the "cloud-like" emitting clumps 
(n, (r) in our notation) and that of the "cloud- like" struc- 
tures' distribution {p*{r) in our notation) on the radial 
distance r from the central black hole. This analysis has 
been applied by the authors mentioned above to the case 
of the well studied and representative Seyfert 1 NGC 5548. 

Combining a photoionization code with an optimiza- 
tion routine. Goad & Koratkar (1998) have found that 
line ratios and variability time scales can be reproduced 
assuming (r) cx Kaspi & Netzer ( |1999| ), through 
the "direct" method of guessing cloud properties and dis- 
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tribution and calculating the resulting emission lines, have 
also found that a density gradient is necessary but their 
best fit is n*(r) oc with 1 < s < 1.5. In our results 
n^: (r) slope is variable in the interval [ri , where star 
envelopes contribute to build up the covering factor. For 
n*(r) profiles shown in Fig.s |[ ^ and ^ a piece- wise power 
law approximation gives logarithmic slopes from 1.5 up 
to 2. These values cannot be used to exclude or to prefer 
one of the assumed mass loss profiles, since they do not 
contradict any of the two previously quoted analyses. 

Recently Korista & Goad ( |2000| ) have tested the "lo- 
cally optimally emitting clouds" (LOG) model, proposed 
by Baldwin ct al. ( |1995D , by comparing the predicted spec- 
trum with that of NGG 5548. One of the results they 
attain is that the power law index of the quantity they 
define as the radial cloud covering fraction must be in 
the range -1.6 to -0.5. In our notation this implies that 
the product R^^^{r)p^{r) must be a function of r which 
is representable by a power law whose exponent should 
be in the range [—1.6, —0.5] mentioned above. Since the 
star distribution function p* (r) that we have adopted has 
a logarithmic slope ranging from -0.7 around r/rg — 10^ 
down to -2 around r/rg = 10^, this comparison favours 
Rext{r) functions flatter than r° ''^. This fact again im- 
plies that "enhanced" and "intermediate" mass loss rates 
are preferable, since in these cases the functional form of 
Rext{r) does indeed fulfill the requirement above (in fact, 
it turns out that it is oc r°-^^ or flatter both for MS and 
RG stars). 



winds; their main result is that these stellar winds can 
be affected in the sense of a mass loss enhancement. Also, 



Both Korista and Goad (200C) and Kaspi and Netzer 
( |1999D derive an estimate for the extension of the BLR. 
For the first ones the maximum BLR extension can be up 
to 200 hght days, while for the second ones it can be up to 
100 light days. In our model, the values for the derived ex- 
ternal BLR radius, r2, for the case of "intermediate" mass 
loss (see Sect. 5.3) are larger than those derived for the 
case of "enhanced" mass loss rates presented in Sect. 5.2; 
however, they are in any case consis tent w ith the esti- 
mate obtained by Korista and Goad (200C). Again, this 
comparison is encouraging for our model, but not definite 
enough to allow for a selection of most appropriate mass 
loss rate values. In any case, the general comparison of 
our model with the global inferred physical conditions for 
the BLR discussed above favours the two configurations 
analyzed in Sects. 5.2 and 5.3, characterized by distance 
dependent stellar mass loss rates and especially by a sub- 
stantial contribution to BLR emission coming from MS 
star envelopes, whose mass loss rate is strongly enhanced 
with respect to their "standard" value. 

Besides the X-ray illumination mechanism that we 
have discussed and examined in Sect. 5.2, other mech- 
anisms have been analyzed in literature, that support 
the physical plausibility of this type of induced and r- 
dependent stellar mass loss rates. Here we just want to 
mention a couple of these works. For example, MacDonald, 



Stanev & Biermann (1991) performed an analysis of the 
effects of neutrino and high energy particle fiux from the 
central source on the stars of the central cluster and their 



Baldwin et al. (1996) have computed the radiative acceler- 
ation on the photospheric layers of stars at BLR distances 
for a number of high-luminosity AGNs; comparing this 
acceleration with the star surface gravity, they conclude 
that all the stars, including main sequence ones, can in- 
deed be affected, thus producing an enhanced mass loss. 
We recall that for the case of r-dependent and enhanced 
mass loss rates that we have discussed in Sect. 5.3, we 
have not explicitly defined the physical mechanism induc- 
ing the enhancement, similarly to what is done in Scoville 
and Norman ( |1988D (see Sect. 5.3), but we just supposed 
one of such mechanisms to be at work. 

Going back to what we have pointed out above, the 
two cases discussed in Sects. 5.2 and 5.3 fall in the same 
"category" of models, characterized by strongly enhanced 
mass loss rates and a significant contribution to the BLR 
emission coming from MS stars. This is interesting from 
two different, but connected, points of view. First, within 
this picture the large number of individual emitting units 
does not imply too large collision rates between the stars. 
In fact, the star collision rate (Eq. (p^) relative to the 
star distribution, p*(r), for which we have computed the 
models whose results are shown in Sect. 5 (the one labeled 
with "1" in Fig. |l]) turns out to be r ~ 1 collisions/yr. The 
total star collision rate is an estimate of the total number 
of stellar collisions that take place in the whole cluster in 
a year; however, to understand to what degree the cluster 
structure is affected by stellar collisions another quantity 
must be introduced, that is the time that it takes for a star 
to be destroyed by collisions with other stars belonging to 
the cluster. Following Begelman & Sikora (1992), we in- 
troduce here the destruction time due to mutual collisions 
for a star in the cluster as 



tcolK?") 



10 



P*{r)h '^{R.Wiry 



(34) 



where the number 10 at numerator is the number of colli- 
sions that a star has to undergo to be finally disrupted; the 
value used above is estimated from the fraction of stellar 



mass lost per collision that Begelman & Sikora ( 1992 ) eval- 
uate as ^ 0.1. In our model the destruction time per star 
turns out to be the same for red giant and main sequence 
stars (since {R*)kg = 10 {R*)ms and /r,g = 0.01 /ms) 
and it varies with the distance r from the central black 
hole. To derive an estimate of its value, we have com- 
puted it in the region of interest for this work, that is at 
a distance which is around a half of the estimated exter- 
nal radius (r2) of the BLR model, i.e. ^ iC^g, obtaining 
tcoi\{x = 10**) ~ 8 X lO^yr. A necessary condition for our 
picture to be consistent is that the MS star evolution time 
is shorter than the star destruction time. Taking into ac- 
count the value of tcow derived above, this condition is 
verified for stars with masses larger than 2.6 Mq which, 
in a star cluster of 3 x 10^ stars, assuming the Salpeter 
initial mass function, constitute ~ 3% of the total star 
number. This percentage of stars can account for the as- 
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sumed quantity of evolved stars in our model. Obviously, 
going closer to the cluster center, icoii decreases and this 
condition can no longer be fulfilled. However, there are a 
few arguments that suggest that the destruction time for 
collisions computed above should be considered as a sort 
of lower limit evaluation. In fact, as Scoville & Norman 
( |1988D argue, ordered stellar motions may result in a lower 
collision rate. In addition, as discussed by these same au- 
thors, stars orbiting on elongated elliptical orbits would 
spend only a short part of their life in the inner region 
of the cluster and this would result in a larger survival 
probability as well. Our present discussion is centered on 
survival conditions for the central stellar cluster as we have 
chosen to model it. Its relation to the AGN lifetime (pos- 
sibly < lO^yr) and to the time at which the nuclear ac- 
tivity switches on with respect to the evolutionary stage 
of the cluster itself would indeed deserve further insight. 
However, in the present context, it is probably just the 
case to mention the possible relevance of this issue and 
we leave its analysis to different work. Back to the orig- 
inal point, even accounting for these last considerations, 
the discussion above shows that the very large values of 
the star number density characterizing the stellar distri- 
bution (r) that we have chosen for our models (the one 
labelled with 1 in Fig.(0)) can be regarded as very close 
to the upper limit for the maintainance of the cluster in- 
tegrity. Indeed, higher stellar densities would result in a 
cluster evolution much faster than what is required when 
the cluster star envelopes are believed to represent the 
structural components of an AGN Broad Line Region. As 
a matter of fact, this is indeed the case when the only stel- 
lar populations forming expanding envelopes suitable for 
Broad Line emission (i.e., contributing to the BLR) are 
those of RGs and SGs. In fact, both RG and SG popula- 
tions amount to only a small fraction of the total stellar 
number in the cluster (i.e., ~ 1%, see Sect. 4.1); as a con- 
sequence, requiring a number density of RG and SG en- 
velopes, as contributing emitting units, sufficiently large 
to justify the observationally inferred covering factors im- 
plies a total number density of stars around 100 times 
larger. This would lead to a destruction time (as estimated 
throug h Eq . (^4|)) unacceptably short, (see Begelman & 
Sikora |l992| ). 

In addition, a second point of interest regarding the 
significant contribution of MS star envelopes to the BLR 
emission, is that, owing to the large number of emitting 
structures, our model can be easily set in accordance with 
the results of the studies of bright AGN spectra, performed 
with cross-correlation techniques, such as those of Arav et 
al. (1997, 1998), that claim to determine a lower limit on 



the number of individual emitting structures of the BLR, 
for models based on discrete emitting units composing the 
region. In particular, for Mk 335, these same authors de- 
rive a lower limit for the number of reprocessing cloud- like 
structures, which is expected to be around 3 x 10®. 

To compare our model to the results inferred by Arav 



Therefore, we have analyzed models for L = 10""^ erg/s, 
that is the luminosity of this source. To model an AGN 
with this luminosity, we have chosen the same central 
black hole mass as that shown in Table 1 of Paper I, 
namely Mbh = 1-12 x 10^ Mq; also, for simplicity, we sup- 
pose that the central star cluster has just the same profile 
p*{r) as the one adopted for the case of luminosity L = 
10'*'' erg/s (see the curve labeled with "1" in Fig. p. With 
these hypotheses, we can test the resulting total number 
of stars within a distance corresponding to the estimated 
external boundary of the BLR model, comparing it with 



Arav et al. (1997, 1998) lower limit evaluation of the num- 
ber of individually emitting units in the source mentioned. 
Our result turns out to be (A^*tot)45 — 3.9 x 10^ stars, 
where we have computed this number using the estimate 
we obtain for (7-2)148 ~ 2.3 x lO^rg, since the total num- 
ber of stars is essentially determined by the number of 
MS stars in the cluster for scenarii of the type described 
in Sects. 5.2 and 5.3. Our evaluation of (A^*tot)45 turns 
out to be w ell above the lower limit found by Arav et al. 
( 1997 , 1998), and this is another encouraging outcome. As 
for other properties of broad line emitting stellar envelopes 
in this same case of "high" luminosity (L — 10'*^ erg/s) 
AGN, these turn out to show behaviours that are simi- 
lar to those we have discussed in Sects. 5.2 and 5.3 for 
the case for L — 10'*'' erg/s. We note in passing that the 
bow shock mechanism for envelope confinement is in gen- 
eral still the most efficient, although in this higher lu- 
minosity case Comptonization confinement (see Sect. 2, 
Eq. (p^)) can be dominant in the inner (r ^ 500rg) por- 
tion of the model broad line emitting zone (that is, closer 
to the central luminosity source, where the radiation flux 
is stronger). 

As a final remark, we note that the fact that the model 
outcome in this scenario can be set in accordance with the 
very large number of discrete broad line emitters in the 
r egion requ ired by the analysis performed by Arav et al. 
( 1997 , 1998) could be of relevance with respect to the BLR 
structure problem, because it would show that discrete 
emitting units models of the BLR a re not necessarily ruled 
out by the analysis of Arav et al. ( 1998 ). 



et al. (1997, 1998) for Mk 335, we have to take into account 



configurations with a more powerful central AGN nucleus. 



7. Conclusions 

In this paper we have analyzed the scenario in which BLR 
emission originates in expanding atmospheres of the stars 
of a dense central cluster embedded in the hot, tenuous 
plasma outflowing from the AGN. Immersing the central 
star cluster in its speciflc AGN environment and investi- 
gating the interaction of the stellar wind with the AGN 
outflow, we consider the shock fronts generated by the 
interaction of the stellar wind and AGN outflow as an- 
other confining mechanism for the stellar envelopes. In 
our model, the nuclear wind structure depends on the fea- 
tures of the dense stellar cluster around the central black 
hole (including the mass loss rates of the stellar compo- 
nents of the cluster itself). At the same time, the physical 
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conditions of the stellar envelopes emitting the broad lines 
depend on the nuclear wind properties. Thus, looking for 
the stellar envelope physical conditions that are appropri- 
ate to reproduce the general characteristics of the BLR as 
they are deduced from observations, we have solved the 
problem consistently. 

We have identified those parameters that significantly 
influence the resulting BLR model properties: these turn 
out to be vq, a stellar wind velocity parameter, the set 
{Mms, M^Q, AIsg} of stellar mass loss rates, and p*{r), 
the star number density distribution. Indeed, we have ex- 
plored different possibilities for building a BLR composed 
of mass losing stellar envelopes, determining, in the very 
specific and extreme AGN central region environment, the 
conditions under which different star types contribute to 
the observed line emission. In summary, we have selected 
two main choices for the BLR structure, namely 1) a 
scenario in which stars are taken as characterized by their 
"standard" (i.e., as known from solar neighbourhood stud- 
ies) mass loss rate, examined in Sect. 5.1, and 2) another 
one in which stars are characterized by an enhanced mass 
loss rate. Both cases, the "enhanced" mass loss rates and 
the "intermediate" mass loss rates, examined in Sects. 5.2 
and 5.3, belong to this second picture. In the first scenario, 
red giant and supergiant envelopes are the foremost con- 
tributors to BLR, while in the second, main sequence stars 
turn out to be an important constituent of the BLR. 

For both scenarios, we have built configurations that 
meet the basic physical conditions for stellar extended at- 
mospheres to be able to contribute to BLR emission. As 
discussed in Sect. 6, the observational requirements are 
much better fulfilled by models built in the second sce- 
nario, allowing for large mass loss rates, with respect to 
those that we can construct starting from the assumption 
of standard stellar mass loss rates. In addition, if we com- 
pare the two scenarios with the further conditions inferred 
from a global analysis of BLR observa tions and with the 
deductions of Arav et al. ( 1997 , 1998 ) on the total num- 
ber of individually emitting structures in the BLR, again 
we find a differentiation in behaviour and outcomes. The 
scenario in which main sequence stars are also endowed 
with winds capable of producing stellar envelopes with a 
non-negligible mass loss rate is therefore preferable. This 
point shows the importance of the choice of the mass loss 
rate in our model. 

Besides the above general considerations, in our model 
the most stringent quantity turns out to be Ctot ■ In partic- 
ular, covering factors that match the range derived from 
observations are not easily obtained. To recover values of 
Ctot as large as those observationally inferred we need to 
suitably choose the star number density and/or the 

ratio M/vq. The star number density distribution of the 
cluster orbiting around the central black hole turns out to 
be crucial in our model. In fact, the observed BLR fea- 
tures can not be recovered for any given p*{r): while the 
required global amount of stars in the cluster is not large, 
the star number density in the central region (i.e. around 



10^^ cm) must be rather high. Also, high mass loss rates 
and low outflow velocities are preferable. 

Even with these limitations, the model nevertheless 
constitutes a consistent structure for a BLR embedded in 
an AGN wind-type plasma outflow, with individual emit- 
ting units characterized by the appropriate values of phys- 
ical quantities for broad line emission. Of course, a more 
detailed comparison of the line shape and intensity pro- 
duced by our model would be necessary to strengthen the 
validity of this approach. As a first step, however, we can 
assert that the results we obtain are encouraging, because 
this analysis has determined strong constraints on the free 
parameters in the model. Such constraints will be useful 
in the derivation of the characteristic line profiles that can 
be predicted by our model. 
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